Integrated guidance and control design by active disturbance rejection method for high-velocity target interceptor with DCS thruster

The present paper proposes a novel integrated guidance and control (IGC) method for engaging with high-speed targets such as ballistic projectiles. considering an extreme short period of terminal engagement due to high relative velocity between target and interceptor, it is particularly important for IGC law to show desirable performance in the presence of various uncertainties (e.g. variation in aerodynamic coefficients) and disturbances (e.g. target maneuver and drag). This article extends the ICG law for mismatched and feedback form equations based on the Active Disturbance Rejection Control (ADRC) method using the back-stepping technique and the Reduced-order Extended State Observer (RESO). The primary consideration is the application of thrusters on the center of mass as the Divert Control System (DCS), along with the daisy-chain technique for control allocation between the fins and thruster commands. Contrary to previous research, the filter and angle measurement error are modeled for the seeker as a crucial parameter to highlight the significance of the thruster. The simulation results indicate the efficiency of the developed method for near-miss or hit-to-kill engagement with tactical ballistic targets. It is shown that the thruster plays a significant role in high-altitude engagements, specifically in the presence of non-ideal seeker. Finally, using the Monte Carlo simulation, it is proved that adding inner loops to the developed technique will not remove the IGC’s advantage over the conventional approach and Non-singular Terminal Sliding Mode (NTSM) guidance law.


Engagement with ballistic targets and its challenges
As a result of the development of ballistic missiles as high-speed vehicles, combating these missiles has been the main challenge for air defense systems.The priority of countermeasures against these targets is interception during the boost phase, followed by exit from the atmosphere at the mid-course phase.This is related to low velocity besides the large radar cross section in an early stage of flight.Later, despite the high velocity, there is no drag or maneuverability due to the lack of atmosphere.Therefore, the trajectory of the projectile will be completely predictable.In these two phases, countermeasures face numerous tactical and technological obstacles such as the need to be in the enemy's area for engagement in the boost phase, or the existence of an extensive and powerful radar network even in other countries.
Consequently, most air defense systems intercept ballistic threats at the reentry phase.During this phase, the target has a high velocity, a small radar cross section, and the ability to change its trajectory.The first two characteristics reduce the interception time, while drag and maneuverability complicate the geometry of the engagement.
To strike ballistic targets, the incident must first be a near-miss or a hit-to-kill one; second, it must occur at high altitudes, where the target cannot maneuver and if one interceptor fails to engage the target, another one still has a chance.Additionally, the cluster or chemical warhead causes less damage to the environment.Loss of aerodynamic efficiency and an increase in the time constant of the interceptor's autopilot are obstacles posed by

IGC and DOBC
In 15,16 ESO estimates target acceleration.Reference 17 focused on ground target with negligible maneuver.It uses ESO combined with back-stepping method.Reference 18 employed a Reduced-order Extended State Observer (RESO) filter and the back-stepping control, demonstrating that RESO has a wider bandwidth than ESO.The actuator rate saturation is not considered in the simulation, and it does not have the ability to engage with high-speed targets with high Zero Effort Miss (ZEM).SMC and super-twisting ESO combination was used for three-dimensional interception, considering the impact angle described in 19 .Due to the type of filter used along with the sliding mode structure, it cannot have the desired efficiency in a short time and at a high speed.To get a faster answer, Non-singular Terminal Sliding Mode control (NTSM) with ESO is used in 20 for intercepting maneuvering target.To avoid chattering, the authors used nonlinear tracking differentiator that complicates the issue.In 21 , RESO combination with the back-stepping and sliding mode framework is also used.The important point of this research is to consider the delay of the actuator.However the desired performance at high speed

Our contributions
As reviewed in Section "IGC and DOBC" most of the research conducted in the field of combining DOBC and IGC has used the structure of back-stepping or SMC due to easy implementation.However what is neglected, is the simultaneous high maneuver target, high closing velocity and high altitude engagement, with multiple disturbances, all required to intercept a ballistic target.These constraints cause the equations to be changed for the use of Divert Control System (DCS) thruster and the requirement of simultaneous commanding to thruster and fins in the final phase arises.Also, the RESO filter is used to estimate the disturbances and uncertainties with a suitable bandwidth and compatible with the cascade control structure.The next issue is that all the cases investigated, use true value for line-of-sight rate without considering the seeker filter and measurement error in the guidance as an important source of error to reach the hit-to-kill interception.In this regard, the main contributions of the present study are summarized as follows: • To overcome the challenges of engagement with a tactical ballistic target, we use an interceptor with a tail and a thruster on the center of mass using the control allocation algorithm in the IGC and DOBC structure.
• Contrary to all existing solutions, seeker dynamics and measurement errors are formulated and implemented.
Due to the short homing phase, this error in the line-of-sight rate poses a significant difficulty in intercepting high-speed targets with a small radar cross-sectional area.• A complete simulation was performed by sweeping ZEMs and different heights to demonstrate the effective- ness of this method versus the conventional method.
It has been mentioned in some references such as 22 that the cascade structure in IGC has destroyed its advantage over the conventional method.In Section "Case 3" this claim is completely rejected and the reasoning is explained.

Organization
The current paper is structured as follows: In Section "Problem formulation" a mathematical model of the thruster is developed for the engagement problem.Section "IGC law design with Reduced order ESO" discusses the design of the controller with an observer.Section "Stability analysis of the closed-loop system" demonstrates the stability of the implemented method, while Section "Seeker filter design" defines seeker dynamics and measurement error.Furthermore, Section "Simulation results" presents the results of detailed comparisons between simulations to illustrate the crucial role of the proposed method.Section "Concluding Remarks" provides conclusion remarks, discussing the effectiveness and difficulties of proposed method for intercepting ballistic targets.

Problem formulation
In this section, the mathematical model of the engagement kinematics is derived.Then the nonlinear dynamic model of the interceptor with thruster is used to develop the integrated guidance and control equations in the pitch plane.As mentioned in the previous section, the IGC system is considered in the homing phase and does not affect other flight phases.Then, the control goal of the paper is described.

Engagement kinematics
The planar geometry of the missile and ballistic target engagement in the inertial system of X I -O I -Z I is shown in Fig. 1, where the missile and target are denoted by M and T, respectively.V t and V m are target and missile www.nature.com/scientificreports/velocities and γ t and γ m are flight path angles, respectively.Also a t and a m are normal accelerations.In addition, R is the relative distance, and is the line-of-sight (LOS) angle.The missile-target relative motion kinematic model is established as 22 : where V m is assumed as constant, i.e., Vm = 0 , and differentiating (1-II ) with respect to time and considering (1-I)-(1-IV ), yield: Remark 1 The term a t cos(γ t − ) is the acceleration perpendicular to the target's line-of-sight.Tactical ballistic targets generally do not have course correction maneuvers or escape maneuvers in the reentry phase, but even these targets can have accelerations due to the presence of fin installation errors 24 .

Remark 2
The term Vt sin(γ t − ) is not considered in most studies in the terminal phase due to low drag, but it should be considered for ballistic targets in equations because of high velocity and high drag.For example, it can be seen in the simulation section that the speed of a tactical ballistic target in the final phase is about 2.7 times that of the interceptor, as a result, its drag force will be about 8 times, and it will be important to consider the drag acceleration even in a short time.However, if the velocity vector of these targets is not in line with the lineof-sight (which is not in most scenarios), a fraction of drag acceleration perpendicular to LOS will be projected, which makes complex near-miss engagement.For this reason, estimating this acceleration and its compensation in the interceptor guidance, significantly affects engagement success.In this study, the observer, estimates both acceleration terms perpendicular to LOS (due to maneuver and drag), which is used in the IGC law.

Nonlinear dynamic model
The nonlinear model of the interceptor with ACS thruster has been presented in 25 .The model of the interceptor with tail fin and DCS thrusters (as shown in Fig. 2) is derived based on this model: where α is the angle of attack, ω y denotes the pitch rate, θ is the pitch angle, m, I YY are the missile mass and pitch moment of inertia, and F z , M y denote the lift force and the pitch moment, respectively.The corresponding expressions are: (1) www.nature.com/scientificreports/where q is the dynamic pressure, S is the aerodynamic reference area, d is the reference length, F th is the thruster force, C z α and C z δ are the lift force derivatives with respect to α and δ .Also δ is the deflection angle for pitch control and C m α , C m δ and C m ωy are the pitch moment derivatives with respect to α , δ and ω y , respectively.Considering Eqs. ( 2) to (4) and defining x 1 = ˙ , x 2 = α , x 3 = ω y , u 1 = F th , u 2 = δ , the integrated model can be achieved as follows: where and Assumption 1 The term C zδ in Eq. ( 4) is neglected because of the low lift of control tails compared to the body lift.
Assumption 2 Both the actuator and thruster have the bounds of (δ, δ) and (F th , F th ) , because of the physical limitation of the actuator and pressure limitation in the gas generator of the thruster.

Assumption 3
In the terminal phase of engagement, Ṙ (by using the seeker data) and R (by using the fusion of seeker and radar data) are provided with acceptable accuracy.Also, the condition of successful engagement is that R is in the range of [0.1, 1] m.

Remark 3
In most references such as 18,22 , the term cos(α) is not considered in Eq. ( 6).This does not seem right because the equation is derived with the assumption of acceleration being perpendicular to the velocity vector, and the acceleration due to C zα and thruster are perpendicular to the body x-axis.
It is observed that the system of Eq. ( 5) is in the feedback form with mismatched uncertainty.These uncertainties are variable with time and functions of state variables.Also, d 1 is the acceleration perpendicular to the target ′ s LOS due to drag or maneuver, while d 2 and d 3 represent the time-varying perturbations caused by vari- ations of aerodynamic parameters and external disturbances i .

Thruster description
Thrusters are typically used in two situations: (a) Attitude Control System (ACS) thruster with less force at a specific distance from the center of mass; This case aims to apply the torque produced by thrusters and rapid rotation of the interceptor to capture the Angle Of Attack (AOA) and, as a result, produce lift force to increase acceleration in the desired direction.In this method, by increasing the altitude, air density decreases and both affect the interceptor's acceleration.(b) DCS thruster with a greater center of mass force; This case aims to generate acceleration in the desired direction for a given period of time.In this class, the thruster acceleration will be independent of altitude.The interceptor with a thruster with four nozzles and the maximum force of T ( 4) on the center of mass can generate an acceleration in a square area by controlling the valves of each nozzle, as shown in Fig. 3.
In this paper, it is assumed that the thruster computer can generate the arbitrary controller force in the pitch plane by changing the valves of its nozzle with good resolution in place, and there is no need for the discretization of the controller's output command.The operation time of a thruster is 1s before the intercept, and it is activated by estimating t go .

IGC law design objective
In this study, the goal is to design an IGC law such that, subject to system (5), a miss distance less than 1 m is achieved against tactical ballistic target.For this purpose, the IGC law should guarantee that the state variable ˙ (x 1 ) is kept close to zero at the end of the terminal phase.Also, one should not be concerned about disturbance d 1 becoming infinity because, as the distance between the interceptor and target decreases to lower than 0.1 m, the simulation is stopped, and destruction occurs in this situation in reality.

IGC law design with reduced order ESO
It is not possible to apply the classical ADRC structure so that all uncertainties are estimated using an ESO filter and then the estimated values are compensated to nullify LOS rate.The existence of mismatched uncertainty and the feedback form structure of the interceptor equation with two inputs, which differs fundamentally from the classical ADRC, are the root causes of this issue.Therefore, a novel concept is employed to deal with these equations by employing the ADRC concept and the framework for dealing with mismatched equations.The observer is created for each category of equations in this technique, which is identical to the ones used in 18,20 .The IGC law in this study is designed using the back-stepping structure and the RESO observer.

Back-stepping based IGC law
A seeker is used to measure the relative parameters of the target and the interceptor.The input of equations defined in the first stage is x 1 d = 0 .Below are three typical modes: 1.The thruster input (u 1 ) is considered zero in this mode and the control structure works to achieve the desired deflection since the thruster activation time corresponding to t go has not yet occurred.In this mode, the first Eq. in (5) changes as follows: The thruster is activated.In this instance, the thruster is initially responsible for providing x c 1 (noting that the thruster time constant is less than fin actuator's one).The deflection command is regarded as zero and the controller structure modifies the thruster force and direction until x 1 approaches the required value if the computed thrust force is less than the maximum thrust value.3. The thruster is activated and the needed thrust force exceeds the thruster's maximum force.In this situation, the angle of attack, which is determined by x c 2 , is responsible for supplying the necessary acceleration difference.The required values to compute the reference command of angular velocity are obtained by sending x c 2 through a differentiation filter.Subsequently, the necessary deflection command is given using x c 3 and ẋc 3 .
All these phases assume that the controller has access to di (the estimation of d i ) instead of d i , which is the out- put of RESO filter with appropriate bandwidth.In the following, the design of the IGC law for the first mode is completely done, and then with the addition of the thruster and use of the daisy chain method, this process is developed for the third mode.The second mode does not require revision because it is only a particular instance of the third mode.www.nature.com/scientificreports/

Step1
The purpose of IGC is zeroing LOS rate ( ˙ ), that is, x c 1 = 0. Equation (8) shows the expected dynamics to meet this demand: where k 1 is the controller coefficient and indicates the convergence speed to the origin.
The first dynamic surface is defined as: Differentiating ( 9) with respect to time provides: where d1 is the estimation of d 1 .The dynamics (8) may then be obtained by determining x c 2 as a virtual command in the form of (11) using the inverse dynamics method.
The objective of the second phase is to choose x c 3 such that the the x 2 state, or the angle of attack ( α ), tracks the command value.This calls for altering the dynamics of the angle of attack, as follows: where Eq. ( 10) is used to calculate α c .A differentiator filter is used to compute αc as follows: where ω is the bandwidth of the derivative filter and "s ′′ is the Laplace notation.Filter ( 13) is used to prevent the explosion of complexity in the analytical calculation.For this computation, there are different approaches, such as using the command filter 26 .Due to the short duration of the final phase and the presence of a thruster, it is not required to employ these methodologies in this study.
By defining the second dynamic surface as: and differentiating it with respect to time, we have: where d2 is the estimation of d 2 .Now, to achieve dynamics in (12), it is suggested to choose the virtual command x c 3 as follows: The controller design without thruster mode is completed after the pitch rate achieves the command value utilizing the suitable elevator.For this reason, the third dynamic surface was defined in the final step as follows: By differentiating ( 17) with respect to time, together with ( 14) yields: where d3 is the estimation of d 3 .Stable dynamics ( 20) is generated to send ω y to ω c y by passing x c 3 through (13) and computing the elevator input as follows: The design created for the first mode is now modified by using the thruster with daisy chain method.The first option in this situation will be to zero the LOS rate by thruster.Since it has a shorter time constant than the electromechanical actuator and does not encounter the rotation rate limit of electromechanical actuators.Then, the difference between the maximum thruster acceleration and the command acceleration is provided by the angle of attack (x 2 ) .Consequently, Eqs. ( 10) and ( 11) are converted to forms ( 21) and ( 22), respectively: www.nature.com/scientificreports/From this stage on, the definition of z 2 is continued as previous.The angle of attack command and therefore the elevator command will be zero if the necessary acceleration is less than the thruster's maximum acceleration.
It is worth noting that in this disturbance rejection-based control structure, α c functions as the guidance loop's output while ω c y acts as the controller's inner loop command.As mentioned in reference 22 , when this control structure is used, the advantage of IGC compared to conventional guidance and control is lost.However, it should be noted in this structure, the dynamics of the interceptor is still considered in the calculation of the guidance command.This differs from the conventional separate designs for the guidance and control loops, which take the dynamics of the interceptor as a point mass in guidance loop.This issue performs better in high-speed engagement geometry changes.Then, the inherent problem in the classical IGC approach should be considered, which ignores the interceptor's inherent longitudinal and angular dynamics difference, as was mentioned in 27 , and led to the development of the partial IGC approach with a structure similar to that used in this study.

Observer design
It was assumed in the previous part that di is provided as the estimator's output.A variety of concepts can be applied while designing the estimator, including ESO, Super twisting ESO, Reduced-order ESO, and high-order nonlinear filters.Each of these approaches has its specific benefit; for example, using super-twisting ESO can ensure the filter's convergence in a finite-time 19 or using high-order nonlinear filters can improve performance when predicting high-frequency disturbances.The RESO estimator, which according to 18 has a greater bandwidth than the ESO, is utilized in this study because of the high speed terminal phase and the lack of high-frequency disturbances.This section reviews how to implement this approach by designing reduced-order ESO as referred in 28 .Since d i is present in the dynamics of all three state variables, an ESO estimator like the classical ADRC cannot be used to estimate the total disturbances; instead, this filter needs to be created independently for each state variable.

Assumption 4
The constant positive values d 1 , d 2 , d 3 exist for uncertainty of the system (5) such that d 1 , d 2 , d 3 meets |d r d i /dt r | ≤ d i ,(i = 1 − 3, r = 0, 1) , i.e. the disturbances and their derivatives are all bounded.
For example, concerning the guidance loop, RESO is as follows: where p 1 , β 1 and d1 represent the filter variable, filter bandwidth, and estimated disturbances in the guidance loop, respectively.In addition, the same procedure is repeated for estimating d2 and d3 , as in Eqs. ( 24) and (25) below: As β increases, the observer bandwidth also increases.Practically, the bandwidth cannot be extended to the intended level due to actuator and sensor data acquisition delays, as well as noise on the sensors.Additionally, the estimation error can be decreased by raising β .Theorem 1 states the relationship between estimation error and β values.
Theorem 1 Considering the observer designed in (23) to (25) and Assumption 4 for system (5), we have: where Proof This theorem has been proved in detail in 18 .It can also be proved using comparison lemma.www.nature.com/scientificreports/By using the comparison lemma, the following inequality is obtained: such that the error of observer is bounded for all t ≥ 0.
Similarly, convergence of other RESO observers (which estimate d2 and d3 ) can be proved.

Remark 4
It should be emphasized that the estimation error will be asymptotically stable by utilizing the observer (23) and adequate adjustment of β if the magnitude of disturbance is constant, i.e. ḋi = 0.

Gain tuning
The bandwidth of the RESO is equal to β .As β increases, the speed of disturbance estimation also increases.Due to the noise of the sensors and also the delay in the system, this value cannot be increased arbitrarily in practice.Controller gains ( k 1 , k 2 and k 3 ), indicate the speed of convergence for LOS rate and tracking for α and q.Also, due to the use of a cascade structure in the control law, the bandwidth of the inner loop must be faster than the middle loop and the outer loop, in order to achieve the desired performance.Here too, the bandwidth of the actuator in practice, makes it impossible to speed up the tracking as much as desired.With these explanations, the following method can be used for the initial gains tuning.First, gain k 3 (which represents the bandwidth of the inner loop) is chosen a little less than the bandwidth of the actuator, and then we choose the bandwidth of the middle and outer loops, respectively, about 2 to 5 times smaller than the previous loop.After that, we should determine the observer's gains.In this regard, the observer gain related to each loop can be selected from 2 to 10 times its control bandwidth, depending on the rate of disturbance changes in that loop.It should be kept in mind that the lack of a sufficient difference between the bandwidth of the inner loop and the actuator or the outer loops can lead to instability or bad tracking.

Stability analysis of the closed-loop system
In this section, the stability of the closed loop system is investigated using Lyapunov theorem.
Theorem 2 Consider the IGC system (5), if Assumption 4 and Theorem 1 are satisfied under the condition that the control gain and observer bandwidth satisfy k i > 0 , β i > 0 , i = 1, 2, 3 , there exists a positive value for k i such that following nonlinear IGC law (29) combined with the RESO estimator can guarantee that the tracking error converges to the origin asymptotically.
Proof In the first step, consider the tracking error as follows: (27) E o is defined by, Matrix P is defined by, By assuming k 1 , k 2 , and k 3 as positive gains which satisfy following inequalities, it can be deduced that matrix P is negative-definite, Keeping in mind, the statement that a matrix is positive (negative) definite if and only if all of its principal minors are positive (negative), the latter is concluded.Now the gains can be tuned as follows: Note that c and r can be set arbitrarily.By using such gains, it can be verified easily that P satisfies following matrix inequality: Hence, V is bounded, and its upper bound is given by, ( 31) V = 1 2

Seeker filter design
Guidance filter plays a crucial role in overall performance of an air defense system 29 .It is well-known that various error sources of onboard seeker, such as low sampling rate besides time-delayed and noisy measurements form the main challenges to achieve a hit-to-kill performance.However, according to authors' knowledge, all existing IGC schemes have assumed an almost ideal guidance filter to utilize true LOS rate value as the measured variable.Keeping this in mind, a two-stage guidance filter is employed in this study to account for exact known engagement kinematics along with an accurate model of seeker error sources.The first filter stage is inspired by 30 which is briefly explained in what follows.
Assume the update rate and measurement delay of seeker to be T s and τ d , respectively.Furthermore, the pointing angle and attributed LOS rate are expressed by ǫ(t) and ω(t) , respectively.The following filter dynam- ics in pitch plane is given: in which r(t) denotes the inertial angular velocity of inner gimbal measured by gimbal's gyroscope.The measurement equation indicates that pointing angle is directly measured by seeker.A classical discrete time filter with gains of L 1 and L 2 is applied to (40) as follows: where and Utilizing the filter dynamics in (41) assures the independence of estimated LOS rates from the interceptor's body angular motions as an important guidance filter performance index 29 .Following straightforward calculations, one can derive the estimated LOS rate dynamic as follows Denoting the standard deviation of pointing angle gaussian noise by σ n , according to (44), the standard deviation of LOS rate estimation is achieved as follows: L 1 and L 2 , are determined in a way which the closed-loop poles of (44) correspond to a standard second order continuous characteristic equation with natural frequency ω c and damping ratio ζ c .It is nice to mention that accounting for range-dependent measurement noise for an active seeker, σ n shall preserve the following equality where σ 0 is the standard deviation of seeker angle measurement error which is a function of signal to noise ratio (SNR) and seeker beam width.Also, R TM 0 is the distance between interceptor and target in the beginning of endgame phase.It is determined according to (46) the decrease of noise level as the relative range goes toward zero.

Remark 5
The main idea behind this section, is to obtain a reasonable model to evaluate the effects of the seeker filter's bandwidth and sensor's measurement noise on proposed guidance law performance, as two significant practical issues.Evidently one can use more sophisticated filter schemes such as those introduced in 29,31,32 to be implemented on seeker's computer.After that, there is a second filter (23) in the interceptor's computer that uses this LOS rate as input to estimate the target's acceleration and uncertainty due to the aerodynamic coefficients.

Simulation results
Numerous numerical simulations were performed, and the results are reported in this section to evaluate the performance of the proposed IGC law in conjunction with the thruster and to demonstrate its ability to intercept high-speed and accelerated targets.This section examines various aspects of this research through four hypothetical situations.All of these scenarios pertain to an interceptor's terminal phase against a tactical ballistic target with a maximum velocity of 2500m/s.The angle between the interceptor's velocity vector and the predicted intercept point (PIP) suggests that the mid-course guidance failed to nullify the zero effort miss (ZEM) at the beginning of the terminal phase.In this circumstance, it was assumed that: The interceptor's specifications are listed in Table 1.
As mentioned above, the interceptor has an active seeker, and the terminal phase range is less than 7km due to issues such as the frequency band of the seeker, its limited power, and the RCS of ballistic targets in that band.Due to the high relative velocity, the homing phase range is highly effective for ZEM ability to compensate.In addition, 0.2 • was chosen as the standard deviation of the angle measurement error (σ 0 ) for this range.
The relative velocity along the LOS and the interceptor-target distance at each instant were utilized to calculate the time to go (t go ) , and the thruster was activated just 1s before the interception.In addition, the thruster had a maximum acceleration of 8 g.
The bandwidth of the differentiator filter is assumed to be 15 rad/s, whereas the bandwidth of the RESO estimators (for all channels) are set as β 1 = 10 , β 2 = 20 and β 3 = 40rad/s.In addition, the controller gains are selected as follows: After passing through the second order transfer function (49), the elevator command was applied to the simulation by passing the rate limit block up to a maximum of 250 • /s and then the saturation block up to a maximum of ±28 • .
where ω n is actuator natural frequency and ζ is damping factor with ω n = 20Hz , ζ = 0.7.
Notably, including the actuator rate limit in the simulation leads to make the case realistic and reduces the controller bandwidth, significantly impacting the engagement outcome.Simulations indicate that actuator rate must be increased for a proper engagement at higher altitudes.Design-wise, this rate should be increased by increasing the interceptor's height and decreasing the hinge moment.
The interceptor was assumed to have a maximum structural load of 22 g due to aerodynamics and thrusters, leading to a restriction of α c .
Finally, miss distance with ZEM calculation is reported at t go = 0.

Remark 6
To calculate ZEM according to reference 33 , the following formula is used: where Ṙ and R are relative velocity vector and relative position vector in inertial frame, respectively.

Case 1
The objective of the first scenario is to engage a ballistic target with the specifications listed in Table 2.This scenario aimed to demonstrate the effectiveness of the proposed IGC method for a dual-controlled interceptor.The terminal phase is begun at an altitude of 12 km and ended at 13 km.The problem carried the following uncertainties and disturbances: 1. Acceleration of the 9g step of the target (as stated previously, the source of this acceleration could be an error in the installation of the ballistic warhead's fins or the presence of a maneuver (47) [ ˙ (0), α(0), ω y (0)] = [ ˙ 0 , 0, 0] www.nature.com/scientificreports/ to change the trajectory of some tactical ballistic missiles); 2.An 8g drag acceleration at the beginning of the scenario; 3. Reduction of aerodynamic coefficients C z α , C m α , and C m δ in the controller relative to simulation equations by 30% , 25% and 25% , respectively; 4. A 7 cm installation error between the thruster nozzle and the interceptor's center of mass.˙ true was supplied directly to ( 10), ( 21) and ( 23) without passing through the filter and adding noise based on the ideal seeker assumption.The kinematics of the engagement is depicted in Figure 4.
The final phase's engagement time was 1.95s, and the obtained miss distance was 0.46 m, as represented in Figure 5.This value indicates the direct interception and success of the scenario.Due to the target's high speed, the velocity vector has rotated approximately 3.5 • due to gravity, maneuver, and drag, which can result in a larger miss distance if the acceleration is not estimated and compensated for in the IGC law.As depicted in Figure 6, the interceptor's velocity vector rotated by 14 • in a short period, causing a significant acceleration to compensate for the error.Before using the thruster, the interceptor utilized all its aerodynamic capabilities to minimize the error.The thruster was then activated 1s before the termination.As the error decreased, the aerodynamic and thruster accelerations decreased too, and in the final moments of the engagement, the interceptor attempted to make the miss distance zero by switching the sign of the acceleration.
Figure 7 depicts the elevator steering so that the LOS rotation rate becomes zero and the pitch rate and AOA track the command in the presence of uncertainties, which is illustrated in Figure 8.
In parallel guidance techniques, such as Proportional Navigation (PN), which are based on zeroing the LOS rate, the guidance bandwidth continuously grows by decreasing the relative range 34 .Consequently, the LOS rate changes drastically at the conclusion of the engagement as depicted in Figure 8. Also, this Figure demonstrates that, given a proper time constant, the AOA and angular velocity can follow their commands.Because the inner loop is quicker, the angular velocity tracking error is less than the AOA.Since tracking differentiator is not used, an initial error with command values is created.
As it is well-established, using a first-order lag and pseudo-differentiator in the scenario is sufficient and does not cause any problems.Before the thruster was activated, the aerodynamic acceleration could not prevent the increase in the LOS rotation rate.Once the thruster was activated, however, this parameter became zero.Finally, Figure 9 depicts both the estimated and actual values of the disturbance in the problem.As it can be seen, the  disturbance leap parameters that occur at 0.85s are the activation of the thruster and the torque produced by its distance from the center of mass.The disturbances caused by the target's drag and acceleration are accurately calculated in d 1 , and the observer's bandwidth is suitable for high speeds and short periods.
In this case, it was demonstrated that the proposed method is suitable for use in destroying a ballistic target.Table 3 presents the miss distance for this scenario in four different modes to emphasize the importance of utilizing the thruster and observer in this engagement.Accordingly, an observer is required to estimate disturbances and uncertainty for the hit-to-kill interception.The impact of the thruster on the interception at high altitudes is another element of significance included in Table 3.Therefore, the absence of a thruster caused a miss distance of 15.3 m at an altitude of 12 km but only 5.5 m at an altitude of 9 km.This is due to a decrease in target velocity resulting from a decrease in interception altitude, required actuator rate, and interceptor time constant.However, it is not recommended to intercept ballistic targets at low altitudes due to factors such as the target's increased maneuverability, the increased risk of destruction, and the presence of cluster warheads.The last thing that is investigated in this case is the effect of the controller gains on the miss distance at the height of 12 km.By changing the gains as k 1 = 0.5 , k 2 = 1 and k 3 = 2 , the simulation results show that the tracking error has increased, but in the presence of the thruster, the miss distance has not changed much and has reached 1.23 m.If the thruster is removed, the miss distance increases to 3.71 m.Also, if the gains are chosen as k 1 = 1 , k 2 = 5 and k 3 = 25 , due to the increase in the bandwidth of the inner loop compared to the actuator, the simulation becomes unstable and the miss distance becomes 136.2 m.The results indicate that a thruster is required to engage high-speed targets, such as ballistic ones when a non-ideal seeker is present.

Case 3
In this case, the Monte Carlo execution for various ZEM's and altitudes is examined to demonstrate the superiority of the proposed approach over the conventional two-loop method.The guidance and control technique described in reference 33 , in which the guidance and control loops are developed independently, serves as a comparison baseline.This method employs the PN law and the three-loop autopilot, as shown in Figure 15.
In this case, the objective was to compare the miss distance between the developed integrated method and the conventional guidance and control one regarding the time lag between the loops.As a result, acceleration was given first to the thruster and then to the three-loop autopilot in the same manner as in the proposed method.In this scenario, a non-ideal seeker was also utilized, as in case 2. No uncertainty was introduced into the problem to reduce the effect of the observer, and the only disturbance was the drag acceleration of the target.Case 3 possessed the same engagement parameters as Case 1.
The three-loop control law is as follows:   where a Z and a z c are interceptor real and command acceleration in pitch plane, k a , k ω , k θ denote the control parameters and design as k a = 2.23 , k ω = 0.26 , k θ = 19.3.
Each point was simulated five times to account for the noise in the seeker measurement, and the mean miss distance at that point was then reported.
As shown in Figs.16 and 17, the performance of the developed integrated guidance and control method (in identical situations like thruster existence and no uncertainty) is superior to that of the conventional method, despite the time lag existence in both methods.It is due to the consideration of the interceptor's dynamics in the calculation of AOA command, and required less time lag in the proposed method.This result validates the claim made in subsection 3.1 As can be seen, larger ZEM's and higher altitudes result in greater increase in miss distance.In addition, the asymmetry of the miss distance around the ZEM is caused by the ballistic target's drag acceleration.

Case 4
In the following explanation, to show the superiority of this method over SMC ones, a comparison with Nonsingular Terminal Sliding Mode (NTSM) guidance law, as the main basis of many new research, has been made.

Figure 4 .
Figure 4.The engagement geometry in Case 1.

Figure 5 .
Figure 5.The changes in path angles, ZEM and relative velocity in case 1.

Figure 6 .
Figure 6.Acceleration of the interceptor due to the angle of attack, thruster activation and the total ones in case 1.

Figure 7 .
Figure 7.The deflection command in case 1.

Figure 8 .Figure 10 .
Figure 8. Changes in the angle of attack, pitch rate and LOS rotation rate in case 1.

Figure 11 .
Figure 11.Changes in the angle of attack, pitch rate and LOS rate in case 2.

Figure 12 . 2 .
Figure 12.Acceleration of the interceptor due to the angle of attack, thruster activation and the total ones in case 2.

Figure 13 .
Figure 13.Changes of dynamic variables in Case 2 without thruster.

Figure 14 .
Figure 14.The deflection command in Case 2 without thruster.

Figure 16 .
Figure 16.The miss distance contour at different altitudes and ZEMs with integrated guidance and control method.

Figure 17 .
Figure 17.The miss distance contour at different altitudes and ZEMs with the conventional guidance and control method.

Table 1 .
Ṙ The interceptor properties in the terminal phase.

Table 2 .
The kinematic properties of the interceptor and the target at the start of the terminal phase.

Table 4 .
Comparison of characteristics of three simulated methods.